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Abstract 

We study the potential to observe CP-violating effects in supersymmetric t\- 
cascade decay chains at the LHC. Asymmetries composed of triple products of 
momenta of the final state particles are sensitive to CP-violating effects. Due to 
large boosts that dilute the asymmetries, these can be difficult to observe. If all 
particle masses in a cascade decay are known, it may be possible to reconstruct all 
momenta in the decay chains on an event-by-event basis even when we have missing 
momentum due to a stable lightest supersymmetric particle. After the reconstruc- 
tion, the non-diluted CP-violating signal can be recovered and gets significantly 
enhanced so that an observation may become feasible. A fully hadronic study has 
been completed to define the areas of the mSUGRA parameter space that may yield 
a 3-0" observation with 500 fb _1 at the LHC. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) is a particularly compelling ex- 
tension of the Standard Model, that may soon be explored at the Large Hadron Collider 
(LHC). It allows one to stabilize the hierarchy between the electroweak (EW) scale and 
the Planck scale and to naturally explain electroweak symmetry breaking (EWSB) by a 
radiative mechanism. The naturalness of the scale of electroweak symmetry breaking and 
the Higgs mass places a rough upper bound on the superpartner masses of several TeV 
and the fits to the electroweak precision data point to a rather light SUSY spectrum [1] . If 
supersymmetry is discovered, many studies will be required to determine the exact details 
of its realisation. One of the interesting issues in this context is CP violation. While the 
observed amount of CP violation in the K and B sectors can be accommodated within 
the SM, another piece of evidence, the baryon asymmetry of the universe, requires a new 
source of CP violation [2-4]. 

The MSSM contains 105 free parameters [5] and a large number of these may have 
non-zero CP-violating phases, see e.g. Ref. [6]. Many of the phases are unphysical in 
the sense that they can be rotated away by a redefinition of the fields. The parameters 
normally chosen to be complex and relevant to this study are the U(l) and SU(3) gaugino 
mass parameters Mi and M 3 , the higgsino mass parameter \x and the trilinear couplings 
of the third generation sfermions Af (f — b,t,r). Hence we have, 

Mi = | Mile** 1 , M 3 = \M 3 \e i(t> \ fi = |//|e^ , A f = \A f \e i(f,A f . (1) 

The two complex parameters that enter the t sector at tree level are A t and fj, and 
in the Xi sector \x and fa. Certain combinations of the CP- violating phases of these 
parameters are constrained by the experimental upper bounds on various electric dipole 
moments (EDMs), see e.g. Ref. [7]. Ignoring possible cancellations, the most severely 
constrained phase is that of fi which contributes to the EDMs at the one-loop level. 
In general for (9(100) GeV supersymmetric masses, |0 M | has to be very small and we 
therefore set M = throughout our study. The phase of A t has weaker constraints as it 
only contributes to the EDMs at the two loop level [8-14]. Here we study the complete 
range of <pA t in order to see the general dependencies exhibited by our observables and 
the luminosity required to observe this within the LHC environment. In principle, <pi can 
also contribute to our observables but in the mSUGRA scenarios discussed in this paper, 
the dependence is weak due to the wino character of the We would like to stress 
that in the chosen scenario experimental bounds from EDMs can be evaded by arranging 
cancellations between various supersymmetric contributions for any value of A t [7,15-17]. 

In general CP phases alter the couplings and masses of SUSY particles, see Ref. [18] 
for a recent review at the LHC. Therefore, in principle we could detect CP-violating 
effects by studying mass spectra, cross sections and branching ratios [19,20]. However, to 
interprete these measurements accurately, we will require high precision and will rely on 
many assumptions of the underlying SUSY breaking mechanism. In addition, all of these 
observables are CP-even and can be faked by a multitude of other parameters. 

In order to make the unambiguous observation of a complex parameter, we need to 
use CP-odd observables. Examples of CP-odd observables include rate asymmetries of 
cross sections and branching ratios. Another possibility, however, are observables that 
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are odd under T-transformations. Applying CPT-invariance, T-odd observables can be 
transferred under certain conditions into CP-odd variables, see Sec. 2.2. These kinds 
of observables can be defined using the triple product correlations of momenta that are 
based on spin correlations of particles, see Refs. [21,22] for a recent review. For the case 
of SUSY at the LHC, we can do this using the final state particles of cascade decays. 

The investigation of triple product correlations within SUSY at the LHC has been 
looked at for various different processes. Stop cascade decays were first studied in Ref. [23] 
and large CP-violating asymmetries were found. The study shown in Ref. [24] was the 
first to specifically examine stop decays in relation to the LHC and significant dilution 
factors were noted when parton distribution functions were introduced. Explicit dilution 
factors and initial estimates of the luminosity required at the LHC were shown for three 
body decays in Ref. [25] and two body decays in Ref. [17]. In addition £ 2 decays were 
investigated in Ref. [26]. b decays have also been looked at in similar studies for two body 
cascade decays [27,28]. 

In Ref. [29] we looked at qg production and decay and studied how to cope with sta- 
tistical limitations and dilution factors in searching for CP-phases in SUSY at the LHC. 
For the present paper, we extend the ideas described in detail in Ref. [29] of momentum 
reconstruction to t production and two-body decays. We further include hadronic, com- 
binatorial and background effects to study whether CP- violation will be observable in the 
t sector at the LHC. 

Regarding a precise measurement of 4>A t at a future linear collider, we are not aware 
of any studies that measure CP-violation in the stop sector directly. However, studies 
have been completed to measure the absolute value of A t and are expected to be accurate 
to within 10% [30,31]. Also, there is the potential to study the CP properties of other 
3rd generation trilinear couplings, namely those of the f sector, (f>A T [32,33]. In addition, 
these studies may also be applicable to the LHC if tau polarisation can be probed in the 
final state [34]. 

For our study, we consider the LHC production process, 

VP Ul- (2) 
Our signal CP-odd observable is generated in the following two body decays, 

ti^X% xl^&N, i^X%, t^b + W. (3) 

where £n and £f denote the near and far leptons respectively. The CP-odd observables 
are built from triple products of final state momentum or reconstructable particles, e.g. 
Pe N ■ (Pt x fry). 

Triple products constructed in this way are not Lorentz invariant but instead depend 
on the intrinsic boost of the produced particle in the laboratory frame. The observed 
asymmetry is maximal when the decay is at rest in the laboratory frame and any boost 
dilutes the observable. Consequently, we decided to use the idea of momentum recon- 
struction to find the momentum of the invisible Xv We are able to perform momentum 
reconstruction for the decay chain shown in Eq. (3) as we have four on-shell mass condi- 
tions which we can solve for the four unknowns of the Xi momentum on an event-by-event 
basis. Once the Xi momentum is known, we can find the rest frame of any particle involved 
in the decay chain and thus measure the maximum CP asymmetry. 
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An important note to make is that the sign of the asymmetry generated by the triple 
product flips if we consider the decay of the charge conjugate t\. Therefore, in addition 
to measuring the triple product we must also determine the charge of the decaying t±. 
Unfortunately we cannot use a leptonically decaying W in this study as we must fully 
measure the t momentum to perform momentum reconstruction. Hence, we rely on the 
opposite ti decay to a single charged lepton final state to tag the charge of both produced 
stops e.g. t\ — ?■ Xii, t — > M~U£. As an aside, charge identification of the process is also 
required to rule out TV-odd observables that can in principle be generated by final state 
interactions at the one-loop level [35], see Sec. 2.2 for more details. We compare the signal 
process with the charge-conjugated decay and if a non-zero asymmetry is observed in the 
combination, it must correspond to a violation of CP symmetry. 

Apart from backgrounds due to hard interactions, measuring asymmetries in a hadronic 
environment is challenging due to the high QCD activity and underlying event that can be 
hard to disentangle from the signal process. However, DO at the Tevatron has succeeded 
in making such a measurement with the like-sign dimuon charge asymmetry [36]. The 
asymmetry is interpreted to originate from the mixing of neutral B mesons and differs by 
3.2 standard deviations from the standard model prediction. Therefore, the measurement 
is a significant hint of a new source of CP-violation. Although this particular measure- 
ment is unlikely to be possible at the LHC due to the pp initial state, it does show that 
if the correct observables are chosen, asymmetry observations are possible at hadron col- 
liders. In addition, CDF at the Tevatron has also recently made the observation of an 
asymmetry in the pair production of top quarks that also hints at new physics [37]. 

We begin in Sec. 2 by describing the process and underlying structure to derive the 
various triple products that can be formed. In Sec. 3 we discuss the momentum recon- 
struction method and its application to the process studied. Sec. 4 gives the analytical 
results of the asymmetries at parton level. Hadron level results are described in Sec. 5 
where we also discuss the effects of standard model and SUSY backgrounds. We also 
find that the method of momentum reconstruction significantly improves the signal to 
background ratio. 

2 Formalism 

2.1 The process studied and the amplitude squared 

At the LHC, the light stop (ti) particles can be produced via pair production, 



which at the LHC will be dominated by the gluon fusion channels, Fig. 1. 

In our study the CP-violating observables are produced in the following decay, 



pp -> htl . 



(4) 



We require the to decay via two, 2-body leptonic channels, 



(5) 



(6) 
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Figure 1: Feynman diagrams for the production process gg — > t\ti. 



where N and F denote the near and far leptons respectively. In addition, we only consider 
events where the t is fully reconstructable and hence decays hadronically, 



t Wb q u q d b . 



(7) 



Using the formalism of Refs. [38,39], the squared amplitude |T| 2 of the full process 
can be factorised into the processes of production gg — > tit* and the subsequent decays 
h — > tx®, X2 ~ * Mn-i £ ~ * Xi^f and t — >■ Wb. We apply the narrow- width approximation 
but include the full spin correlations for the production and the decay of the intermediate 
particles, t\, x®, £ an d t. The use of the narrow-width approximation is appropriate since 
the widths of the respective particles are much smaller than the masses in all cases. The 
squared amplitude can then be expressed in the form, 

|T| 2 = ^ACfOI^A^I^A^I^WI^f^jIp^D^D^D^ 

a=l I 

+ E xfm)xi(xi)xut)D(£) 



b=i 



(8) 



a,b=l 



where a, b = 1, 2, 3 refers to the polarisation states of the neutralino x° and top quark t. 
In addition, 

• A(tx), A^), A(£) and A(i) are the pseudo-propagators of the intermediate parti- 
cles which lead to the factors E^/m^T^, E^o/m^oT^o, E^/m^T^ and Et/m t T t in 
the narrow-width approximation. 
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• P(titi), P{ty$), D {■)$), D(£) and D(t) (Appendix D) are the terms in the production 
and decay that are independent of the spin of the decaying neutralino and top, 
whereas, 

• £p(x°), £$,(<), £p(x°0 and E b D (t) (Appendix D) are the spin-dependent 
terms giving the correlations between production and decay of the x\ an d t. We 
follow the formalism and conventions described in Ref. [39]. 

• It must be noted that the slepton £ produces no spin correlation term in the ampli- 
tude since it is a scalar. 

Explicit expressions are given in Appendix D. 
2.2 Structure of the T-odd asymmetry 

As shown in the CPT-theorem [40,41], relativistic quantum field theories with usual spin- 
statistics relations have to be invariant under a CPT-transformation. This invariance 
guarantees that the masses and also the total widths of particles and antiparticles are 
the same. Since a true T-transformation is anti-unitary, which exchanges the initial 
and the final states, it is useful to study 'naive' TV-transformations for collider-based 
experiments. The definition of Tjv-transformations is to apply T-transformations to the 
initial and final states but without interchanging them. The unitarity of the S-matrix 
leads in the absence of re-scattering effects (i.e. in leading order in perturbation theory, 
no final state interactions (FSI) and no width effects) to a conservation of the scattering 
amplitude under a CPTjv-transformation [35]. 

It is therefore useful to categorise CP-violating observables into T^-odd and Tjy- 
even observables. CPTjv invariance implies that a Tjv-odd observable is also CP-odd in 
the absence of re-scattering effects. However, in case re-scattering effects contribute, i.e. 
CPT^v 7^ CPT- invariance, a T^r-odd signal may be caused by such re-scattering effects 
and does not necessarily imply CP-violation. 

For all our observables we require that we know the charge of the decaying t\ and can 
therefore distinguish the particle and anti-particle. Hence we can combine the process 
with the charge-conjugated decay to make an unambiguous observation of CP- violation 
via T^r-odd observables. 

In general, it is therefore important to classify all terms of the corresponding amplitude 
squared, eq.(8), with respect to their Tjy-odd or T^-even character. Only the products 
that contain a Tjv-odd contribution will lead to CP-odd violating observables: 

• The spin-independent terms introduced in the previous section, P(ti£i), Pftx®), 
D(xX)-> D(£), D(t) do not cause any Tjy-odd terms. 

• The spin- dependent terms, X P (x°), £p(t), s p s i)(X2)> 4W» however, often 
can be divided up into Tjy-even and Tjy-odd terms, depending on the processes 
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studied. In our case, a sequence of 2-body decays, we can only split Sp(x^) = 
^P^evenix^t) + ^p^odJo&t) an d other spin-dependent terms only lead to T/v-even 



terms 1 . 

• Therefore, the TV-odd term in the amplitude is, Y? a ,b=i Zf,oM(xit)Z a D (x%)Z b D (t)D(£). 

When we contract the spin indices of the t and x° an d evaluate the TV-odd contribu- 
tion, we find that the following covariant product appears in the amplitude, 

ZtM°2t)Zl(xi)Z b D (t) ~ ie^s^(xl) P ^/' p m x (Pe N s a )( P[b , w] s b ), (9) 

~ K»vpaP\$t N PwP C t , ( 10 ) 

where T# t0dd , ^oixl) and £$,(*) are given by Eq. (69), Eq. (73) and Eq. (75), respectively. 

The above equation is multiplied by the imaginary part of the coupling, Eq. (71) that 
contains terms from both the t, Eq. (39), and x°, Eq. (50), mixing matrices. Hence, any 
complex phases contained in those mixing matrices will yield CP-violating effects that 
can be seen in an observable that exploits the covariant product. We can now expand 
the Lorentz invariant covariant product in terms of the explicit energy and momentum 
components, 

£^P«P U x° 2 Pi N PwPt = E x% Pzt ■ (Pw x pt) + E w pt ■ (p~| x p^) (11) 
-E £n Pw ■ (Pt x p$>) - E t p~| ■ (p^ x p~w) . 

The first term in Eq. (11) shows the CP sensitive triple product that can be measured 
from final state momenta. However, this triple product is not Lorentz invariant and 
consequently can vary in both magnitude and sign in different reference frames. If we are 
in the rest frame of the x\ though, 

£nup*P^oPe N PwPt — > m x% ■ (Pw x pt) , (12) 

the resulting asymmetry, Eq. (15), is uniquely defined since all other terms of the covariant 
product vanish as p^o — > 0. 

Hence we see that triple products of momenta, can be used as Tjv-odd observables. 
In this paper we find that the triple products most useful to study are, 

Te N = Pi N ■ (pw x p t ) , (13) 
%t = Pb- (Pe+ xpe-)- (14) 

where £ + and £~ are the two leptons produced in the cascade decay. For the triple 
product, Eq. (14), the identification of near and far leptons is not required as is explained 
at the end of this section. 

The T-odd asymmetry is then defined as, 

N T+ - N T _ /sign{7»|r| 2 dUps 

N T+ +N T _ J|T| 2 rflips ' 1 ' 



1 This is different if 3-body decays are studied, see Ref. [25]. In that case spin-dependent terms from 
both the production Sp (x^i) as well as from the 3-body decay Y.^xV) l ea d to CP-odd contributions. 
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where / = £n or ££, rflips denotes Lorentz invariant phase space and Nj- + (Nj-_) are 
the numbers of events for which T is positive (negative). The denominator in Eq. (15), 
J* |T| 2 (i lips, is equal to the total cross section. 
We then define, 

Ae N = A t {T £n ), A u = A T (T U ), (16) 

where Ae N is the asymmetry from the triple product % N and Au is the asymmetry from 
the triple product Tit- 

As stated above, whilst the covariant product is Lorentz invariant, the triple products 
are not. However, we can see that for the triple product in Eq.(13), the rest frame of the 
x\ and the i-y are equivalent since (pf = p^o + p t ), 

(17) 

For the triple product 7«, Eq. (14), the covariant product can be re-expressed in the 
following form (exploiting momentum conservation, p^o = Pi + pe N , Pi = pt F + p^o, 
Pw =Pt+Pb), 

^upaP^oPe N P P w Pt = (pe F + P^P\ N PwV U b • (18) 

We now see that we have effectively two covariant products, one which contains the 
momentum of the Xv I n general, triple products containing the momentum of the far 
lepton will be lower as the far lepton is not directly correlated with the spin of Xt 
Nevertheless, we can exploit and maximise the triple products originating from Eq. (17) 
and Eq. (18), if we know the momentum of the unstable particles in the decay chain. This 
can be provided by the momentum reconstruction procedure described in the following 
section. 

Changing the decaying t\ to a t\ or changing the charge of the near lepton £ N reverses 
the sign of the covariant product. Consequently we have to know the charge of both the 
t\ and the £n, otherwise any asymmetry will cancel. The charge of the t\ can be found 
by demanding that the opposite cascade produces a single lepton and thus a tri-lepton 
final state. We distinguish the near and far leptons using the momentum reconstruction 
technique, Sec. 3. However if for some reason the leptons cannot be identified, we can 
still use the triple product 7~a, Eq. (14). No lepton distinction is required as exchanging 
the near and far leptons has an extra sign change that cancels the change produced by 
the charge exchange. 



3 Momentum reconstruction 
3.1 Dilution effects 

The triple product that is constructed from momenta in the laboratory frame suffers from 
dilution factors (~ 4) at the LHC. This is due to the lab frame being boosted with respect 
to the rest frame of the X2 or see Eq. (17), for a more detailed discussion see Ref. [25]. 
It results in a considerable reduction in the maximum asymmetry observable when we 
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Figure 2: The asymmetry At, Eq. (15), as a function of the stop momentum, \p^\, in the 
laboratory frame. The solid line is the asymmetry for the triple product % N , Eq. (13), 
and the dotted line is for the triple product %e, Eq. (14). The respective masses are given 
in Tab. 2, Tab. 3 and Tab. 4. 

introduce the parton density functions (PDFs) which causes an undetermined boost to 
the system. Fig. 2 shows how the asymmetry is diluted in the laboratory frame when 
we produce the t\ with varying initial momenta. If we were able to reconstruct the 
momentum of the t±, we could perform a Lorentz transformation of all the momenta in 
the triple product into the t± rest frame and potentially recover the full asymmetry. 

3.2 Reconstruction procedure 

We are able to reconstruct the xl f° ur momentum by reconstructing the following two 
body decay chain in full, 



Assuming that all the masses in the decay chain are known, the kinematics can be 
fully reconstructed using the set of invariant mass conditions, 



where p denote the four momenta of the respective particles. 

We see that with the four equations we have enough information to solve the system 
and find each component of the Xi f° ur momentum. A solution to the above set of 



t-tt + xl-tt + P + e^-n + xl + ^ + et- 



(19) 




(Pl± +Pi%f = (Px l l + Pi± +^) 2 , 
(Pxl + Ptf = i.Px\ + P4 + Pe* + Ptfi 



(20) 
(21) 
(22) 
(23) 
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Figure 3: The process studied for momen- 
tum reconstruction. 



equations is presented in Ref. [42] and we outline the procedure here. We first expand 

(IT f± 



the Xi momentum in terms of the final state momentum of the tp, £^ T and t, 



(24) 



In order to derive a system of 3 linear equations for the unknowns a — c, we calculate 
Px° 'Vi F i Px° ' P £ n an d Px° ' Pt- Inserting Eq. (24) and exploiting Eqs. (21-23) we form the 
system of equations, 



M 




( 



\{mlo 
i(mf-m|o) 

2 ™2 \ 



mf) + EjoEi 
+- Pe F ■ Pit, 



E x oE eN 



(25) 



mr 



in 



+ Pi F ■ Pt + Pe N ■ Pt + E x oE t 



where, 



M 



~te F ■ ~fe F 
~fi N ■ ~te F 
~ft ■ ~te F 



~fe F ■ ~fe N 
Pt ■ Pe N 




-ft^ 



(26) 



We invert the matrix M. to find solutions for a, b and c in terms of constants and E^o. 
The on shell mass condition for the x?> Eq. (20), can then be expressed as, 

(a, b, c)M 




(27) 



We solve the above quadratic equation, to find two solutions for E x o. These solutions are 
then substituted back into Eq. (24) to find all components of the t\ momentum on an 
event-by-event basis. 
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3.3 Challenges from multiple solutions 

We encounter a complication in the reconstruction as Eq. (20) is quadratic in p^o. Conse- 
quently we have two solutions for p^o , for each reconstructed event but we have no extra 
information in the single decay chain to determine which solution is physically correct. 
As we cannot distinguish which of these solutions corresponds to the physically correct 
configuration, we need to analyse both. Therefore, we calculate the t\ momentum for 
both configurations and boost all final state particles in the triple product into the re- 
constructed ti rest frame. If the sign of both triple products are the same then the event 
is recorded but if the sign of the triple products are different, we discard the event since 
we cannot know which of the reconstructed solutions is correct. The method has the 
disadvantage that we lose events and therefore statistical significance. However, we find 
that the asymmetry can actually rise 1.5%) as most of the events removed have small 
triple products and events with a small triple product lead to smaller asymmetries. 

The procedure is essentially a cut designed for the triple product correlation observ- 
ables. Events with an ambiguous triple product sign will significantly dilute the asym- 
metry and reduce the statistical significance of any CP-violating observation. Therefore, 
they must be removed from the sample. The disadvantage of the cut is that it makes an 
actual measurement of the CP-violating phase more involved. A comparison would have 
to be performed between a Monte Carlo simulation and the real data for a measurement to 
take place and may induce new errors. However, we believe that an actual determination 
of the phases at the LHC will be challenging and the method presented is more designed 
to establish the presence of CP-violation in SUSY. 

When performing the momentum reconstruction at the LHC we have additional prob- 
lems from multiple solutions that come from combinatorial effects in the event. First, to 
complete the reconstruction we need to correctly identify the near and far lepton in the 
decay chain Eq. (19), if we wish to compute the triple product Te N , Eq. (13), although this 
information is not required for the triple product Tu, Eq. (14). We find that in 20% 
of events the wrong assignment of near and far leptons satisfies the kinematic equations 
Eq. (20)-(23) and produces two extra solutions for the momentum of the Xi m addition to 
the solutions found from the correct configuration. In addition, we always require a third 
lepton in the event coming from the opposite decay chain to correctly identify the stop 
charge. For example the lepton produced in the decay chain i\ — > x~b, xT — > t~ + X, 
where X are other neutral decay products . If this lepton is of the same flavour as those in 
the triple product decay chain there is a small chance that it can also reconstruct the 
All of these combinatorial issues are removed by again demanding that all calculated triple 
products are of the same sign and discarding any events where opposite sign solutions 
occur. 

Further combinatorial issues occur with the reconstructed top in the event. Firstly a 
second b is always present in the opposite decay chain and this can occasionally combine 
with a reconstructed W to give a fake t. The opposite decay chain also can contain 
extra quarks that can produce more reconstructed tfs. Finally, the parton shower can 
sometimes radiate hard gluons that are also seen as extra jets and further complicate 
the combinatorial problem. Whenever extra t quarks are found that satisfy the event 
kinematics, we perform the same procedure as for combinatorial leptons. Triple products 
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are calculated for all reconstructed rest frames and only events, that yield the same sign 
for all the reconstructed triple products, are recorded. 

3.4 Mass measurements 

As mentioned above, we assume that the masses of all the SUSY particles in the decay 
chain will be known. However, for the majority of our equations in Eq. (25), we actually 
require the difference between various m 2 's in the decay chains and not the absolute mass. 
At the LHC, the established way of measuring the SUSY spectrum is via mass end-points 
(see Ref. [43] and references therein) and this method will measure these mass differences 
with high accuracy (9(1%). 

The on-shell mass condition for the Xi requires the absolute mass scale and this should 
be measured at the LHC to a precision of better than 10% [43], for low mass scenarios 
similar to the phenomenology presented in this paper. As an extra check on the numer- 
ical stability of the reconstruction procedure, up to 20 GeV absolute mass errors were 
tested on the absolute mass scale of the decay chain as a conservative estimate. This had 
a negligible effect on the reconstruction efficiency and the CP-asymmetry and is there- 
fore not considered to be a problem. In addition new methods have been proposed for 
measuring the sparticle masses from the kinematic invariants directly [42,44-48]. These 
methods also use the mass invariants on an event-by-event basis but use this information 
to reconstruct the masses of the particles in the decay chain. Therefore, these methods 
are directly measuring the inputs we require for Eq. (20) and Eq. (25). We then use the 
output from these methods to reconstruct the momentum of the Xi on an event-by-event 
basis. Reviews of all the major mass reconstruction methods proposed for the LHC are 
given in Refs. [49,50]. 

4 Parton level results 

In this section we analyse numerically the CP-asymmetry at the parton level, with the 
inclusion of parton distribution functions, while in Sec. 5 we complete a hadronic level 
study to estimate the effect in a realistic environment and the discovery potential at 
the LHC. In particular, we focus on a specific mSUGRA parameter point, Tab. 1, at the 
parton level before discussing more general low mass mSUGRA scenarios for our hadronic 
study. 

4.1 Chosen scenario: spectrum and decay modes 

We choose for this study the mSUGRA scenario shown in Tab. 1 with an added CP- 
phase to the trilinear coupling 0^ t . Although the value of the trilinear coupling is zero 
at the unification scale in this scenario, the renormalisation group equations generate a 
value of, A t = —391 GeV, at the weak scale. The spectrum at the electroweak scale has 
been derived using the RGE code SPheno 2.2.3 [51] and the masses of the gauginos and 
scalars are shown in Tab. 2, Tab. 3 and Tab. 4 respectively. Using the low energy soft 
SUSY breaking parameters and the phase of the trilinear coupling 0^ t , we calculate the 
masses and mixing of the tj's, see Appendix A for details. 
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Parameter 


m Q 


mi/2 


tan (3 


sign(^) 


A 


Value 


65 


210 


5 


+ 






Table 1: mSUGRA benchmark scenario (masses in GeV). 



Particle 






A3 


A4 








Mass(GeV) 


77.7 


142.4 


305.1 


330.3 


140.7 


329.9 


514.1 



Table 2: Masses (in GeV) of the gauginos calculated by SPheno 2.2.3 [51]. 



Particle 


mil 


m h 


m h 


m h 


m <j dL 


m qdR 


™^ 




Mass(GeV) 


345.7 


497.8 


443.4 


466.0 


484.7 


465.2 


478.7 


464.9 



Table 3: Masses (in GeV) of the squarks calculated by SPheno 2.2.3 [51] except for the 
ti which were calculated at tree level for the phase <pA t = |§7r|. 



Particle 






m f2 




Mass(GeV) 


163.4 


110.8 


164.9 


108.0 



Table 4: Masses (in GeV) of the SUSY sleptons calculated by SPheno 2.2.3 [51]. 



For the presented analysis to work, we require the SUSY spectrum to have the following 
mass hierarchy, 

mi x — m t > m^Q > m^± > m^,o, (28) 

to allow for full momentum reconstruction. This hierarchy is often a feature in the 
mSUGRA parameter space. In addition we concentrate on scenarios with a light stop 
as the study is statistically limited and consequently we examine cases with a large pro- 
duction cross section. 2 

The feasibility of the method at the LHC depends heavily on the integrated luminosity. 
For this reason we look closely at the predicted cross section of the asymmetry decay chain, 

a = a(pp tJl) x BR{h i$) X BR(x° 2 I^) x BR(P y&P) x BR(t q u q d b), 

(29) 

2 Since this paper was submitted, the particular parameter point, Tab. 1, has been excluded [52,53]. 
However, the exclusion is derived from the gluino, first and second generation squark masses but the stop 
masses have far weaker bounds [54]. Therefore, if we do not restrict ourselves to a mSUGRA parameter 
space, the study is still valid. 
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Parameter 


Value 


BR(h -)• X°it) 


34.6 


BRfc xlt) 


7.5 


BR{h -> x+fc) 


50.1 


Bi2(f~i 


7.8 




11.6 


BR(xt -> r+z/ r ) 


95.1 


a(pp -)• Ml) [pb] 


3.44 



Table 5: Nominal values of the branching ratios (in %) for various decays calculated in 
Herwig++ [55,56] with phase 4>A t = ||vr|. In the last row, the calculated cross section for 
stop pair production at the LHC with y/s = 14 TeV at leading order (LO) from Herwig++. 



and the relevant values for our scenario are shown in Tab. 5. In our study we also need 
to identify the charge of the t\ in the opposite decay chain and this is possible when the 
decay products contain a single lepton (any number of jets are allowed). We see that 
the dominant production of single leptons from t\ decays is via the channel t\ — > xtb. 
However, as only the right sleptons and the bino-like Xi are lighter than the wino-like 
Xi, the decay of the xf is via mixing terms or Yukawa couplings and hence the decay 
BR(xt — > t~i^t) dominates, Tab. 5. For this reason we find that our study is far more 
promising if r identification is possible. We compare results where r identification (with 
a 40% efficiency in the hadronic channels) has and has not been used in Sec. 5. 

4.2 CP asymmetry at the parton level 

We start by discussing the dependence of the parton level asymmetry on <pA t i Eq. (15), 
for both the triple products % N and Tee, Eqs. (13), (14). In order to see the maximum 
dependence upon <f>A t , we reconstruct the t± at rest and calculate the triple product in 
this frame. It should be noted that the asymmetry is obviously a CP-odd quantity, Fig. 4. 

We see from Fig. 4(a) that the largest asymmetry occurs for the triple product % N) 
which attains l^^lmax ~ 15% when <pA t ~ 0.87T. For the triple product 7~u, the asymmetry 
is smaller, |^4«| max ~ 6.5%, because the 'true' CP triple product correlation is only 
partially measured, see Sec. 2.2. 

If we now include the dominant production process at the LHC (gg — > hti) and 
relevant parton distribution functions (MRST 2004L0 [57]), we see that the asymmetries 
are significantly diluted, Fig. 4(b). The asymmetry for the triple product 7~e N , drops from 
l^jjmax ~ 15% to l^jjmax ~ 4.5% and the reduction is due to the boosted frame of 
the produced ti as discussed in Sec. 2.2. For the triple product Tee, the reduction in the 
asymmetry is far less, from |^4«| max ~ 6.5% to |^4«| max ~ 3.8%. This is because the triple 
product, relies on the £f being correlated with the £ by the intrinsic boost of the x-L ^ 
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(b) 



< 



<t>A t /^ 
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y/s = 14TeV 
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Figure 4: (a) The asymmetry At, Eq. (15), in the rest frame of t\ as a function of 0^ t . 
(b) The asymmetry .At, Eq. (15), in the laboratory frame as a function of <pA t & t the LHC 
at 14TeV. The solid line is the asymmetry for the triple product % N , Eq. (13) and the 
dotted line is for the triple product Tee, Eq. (14). 



system which already has a boost, even when the t\ is at rest. As the t\ becomes boosted, 
the boost of the I system becomes proportionally less so, as the momentum of the 
t\ is distributed throughout the decay chain. The difference in the dilution of the two 
asymmetries with t momentum can be seen in Fig. 2. 



5 Hadron level results 

In order to estimate the potential for observing CP-violating effects in t\ decays at the 
LHC more realistically, we perform the analysis at the hadronic level. We use the 
Herwig++ [55, 56] event generator to calculate all the matrix elements in the process, 
the initial hard interaction, the subsequent SUSY particle decays, the parton shower and 
the hadronisation. An important feature of Herwig++ is that it calculates the spin corre- 
lations [58] in the SUSY cascade decay and allows the input of complex mixing matrices. 
Consequently, the triple product CP-asymmetry can be automatically calculated within 
Herwig++. 

We also include both Standard Model and SUSY backgrounds in the analysis to un- 
derstand how well t\t\ production can be isolated at the LHC. We find that after applying 
basic signal identification cuts and the more complicated momentum reconstruction, virtu- 
ally no Standard Model background contributes, Sec. 5.2. However, the SUSY background 
presents more of a challenge and new cuts have to be introduced to improve the signal 
to background ratio, Sec. 5.4. Even after cuts, the SUSY background can remain prob- 
lematic but if the dominant contributions are known, the backgrounds can be partially 
subtracted. 
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5.1 Cuts used and signal identification 

The hadronic analysis of the produced events has been performed within the program 
Rivet [59,60]. We used the anti-k t [61,62] jet algorithm with R=0.5 and applied the 
following acceptance cuts, 

• pn. > lOGeV, 

• p Tji > 20GeV, 

• invariant mass of opposite sign same flavour (OSSF) leptons: M e + e - > lOGeV, 

• |%| < 2.5, 

• \Vji\ < 3-5, 

• lepton jet isolation, AR = 0.5, 

• 6-tag efficiency = 60% [63], 

• hadronic r-tag efficiency (whenever used) = 40% [63]. 

To identify the events we demand three charged leptons in the final state, so that we 
can correctly identify the charge of each t\ produced in the event, c.f. Sec. 2.2. In addition, 
we demand that a pair of these leptons are OSSF as is the case for light leptons from \\ 
decay. Whenever a t\ decays in our scenarios a b is produced and therefore we require 
at least one 6-tag in the final state (in principle we could require 2 6-tags including the 
opposite decay chain but we lose 40% of events due to b-tagging efficiency). On top of the 
b we require at least 2 more jets to be found in the final state so the full reconstruction 
of the t is possible. As all of our triple products and reconstructed momenta need a t, we 
require at least one hadronic t to be reconstructed. For this procedure, we first demand 
that 2 jets (not 6's) reconstruct a W ± (70GeV < Mjj < 90GeV). We then impose that a 
reconstructed W ± and one b jet reconstruct a t (150GeV < M w ± b < 190GeV). 

Once these cuts have been passed we then perform the kinematical reconstruction 
shown in Sec. 3 with any t's and OSSF leptons found in the final state. If the particles 
satisfy the kinematic constraints Eqs. (20)- (23), we will have at least two different solutions 
on an event-by-event basis for the momentum of the Xv F° r each solution, the relevant 
rest frame triple product is calculated and only if all the signs of the triple products agree 
the event is accepted. 

5.2 Standard model background 

The following standard model backgrounds were produced with Herwig++: it, Drell- 
Yan (via 7 and Z), W+jet, Z+jet, WW, WZ, ZZ and Wj. In addition, we generate 
ti£ + £~ events with MadGraph [64] and then use Herwig++ to perform the parton shower 
and hadronisation. We find that after we produce an equivalent luminosity of 500 fb -1 , 
the only background to pass the event selection is ti£ + £~ with the very low rate of 0.03 
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0.5 1 1.5 2 



0.5 1 1.5 2 



Figure 5: (a) The asymmetry Ag N , Eq. (15), for the decay chain shown in Eq. (30)- 
(31) as a function of <p^ t at the hadronic level after momentum reconstruction has been 
performed, (b) The branching ratios: ti — > xfb (black solid), t± —> xt^ ( re d dotted), 
ti X\t (purple dashed), t\ — > \\t (blue dash-dot). 

events/ib -1 after kinematical reconstruction. This corresponds to only « 1% of the signal 
process for our particular scenario. 

Although the above result is encouraging, it must be stated that our analysis contains 
no jets mis-identified as leptons. As the dominant standard model contributions produced 
by Herwig++ only contain two hard leptons in the initial process, the lack of a trilepton 
signal is not surprising. However, we do not expect major problems from standard model 
backgrounds if we limit the study to leptons from the first and second generation, ti can 
be expected to provide the largest background when both W ± decay leptonically and an 
extra lepton is produced from a h or a mis-identified jet. Even when this occurs though, 
we still require an additional two hard jets in the event that have to combine with a 
b to form a t. Moreover, the final state then has to fulfil the reconstructed particular 
kinematics of our signal and finally all the calculated triple products have to agree. 

To improve the statistical significance of our analysis, we also investigated the possi- 
bility of using r-tagging in the opposite decay chain to that of our signal. In this analysis, 
we now change the original trilepton signal to a first or second lepton OSSF and addi- 
tional hadronic r. The mis-identification of a jet for a r is much higher than for the other 
leptons and the standard model backgrounds may now become an issue [63]. However, 
this analysis is postponed to future studies. 

5.3 Stop Production 

We begin by studying t\t\ production along with the following decay chain, 



*i -> X 2 t -> Xi e+e 3u3d b , 
i*i Xit -> Xiff^b. 



(30) 
(31) 
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tit\ 


9 A 


Herwig++ LO (pb _1 ) 


3.44 


75.8 


Prospino LO (pb _1 ) 




76 7 +24 - 8 

'°- ' -17.3 


Prospino NLO (pb _1 ) 


5-04^;^ 





Table 6: Cross section at the LHC with y/s = 14 TeV production channel tit\ and coloured 
SUSY production for both leading order (LO) and next-to- leading order (NLO). All cross 
sections were calculated using Herwig++ [55,56] or Prospino [66-68]. The errors indicated 
next to the Prospino cross sections relate to varying the factorisation and renormalisation 
scales from 0.5mr — > 2m,- . 



to test the momentum reconstruction procedure. The above decay chain is the cleanest 
signal process from a combinatorial point of view. We find a reconstruction efficiency 
of ~ 5% for this particular topology after cuts and the requirement for same sign triple 
products. The decay chain Eq. (31) has a single lepton in the final state allowing us to 
tag the charge of both the ti and t* in the process. 

For the CP-asymmetry we now concentrate purely on the triple product % N , Eq. (13), 
calculated in the reconstructed rest frame of the ti, as this is the observable with high 
significance at the LHC. Fig. 5(a) shows that there is virtually no dilution when we move 
to the hadronic level and the maximal asymmetry stays at l^jjmax ~ 15%. In fact, 
the hadronic level reconstruction does induce a degree of dilution, 1.5% but this is 
cancelled by our procedure of removing opposite sign triple products which enhances the 
asymmetry by a similar amount, c.f. Sec. 3.3. 

In order to estimate whether it is possible to observe a CP-asymmetry in t\ decays 
at the LHC, we need to calculate the statistical significance of any result. We assume 
that Nj- + (Nj-_), the numbers of events where T is positive (negative) as in Eq. (15), are 
binomially distributed, giving the following statistical error [65], 

A(^ T ) stat = 2a/ e(l — e)/N , (32) 

where e = N r+ /(N T+ + N T _) = \(1 + At), and N = N T+ + N T _ is the total number of 
events. Eq.(32) can be rearranged to give the required number of events for a desired 
significance. 

The total cross section used to calculate the statistical significance of any result in 
this paper has been calculated using Herwig++ at the leading order (LO) for consistency. 
However, next-to-leading order production cross sections are available and have been 
calculated using Prospino [66-68], cf. Tab. 6. We see that in general the cross sections 
at NLO are higher than those at LO suggesting that the effective luminosity at the 
LHC will be more optimistic than those shown in the following results. In addition, the 
factorisation and renormalisation scale uncertainties are shown that indicate an estimate 
of the underlying theoretical uncertainty. 

Due to the phase dependence of both the t\ branching ratios, see Fig. 5(b), and 
production cross section, the statistical significance for different values of <f>A t cannot be 
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Figure 6: Pure tit\ production, all decay channels included, see Tab. 5 for branching 
ratios for the specific parameter point and Fig. 5 for how these alter with <pA t ■ t tagging 
is included in both plots, (a) Asymmetry, Ae N , at reference point with 3a-luminosity 
lines shown, (b) Minimum luminosity required for 3<T-discovery in tan/3, A plane (at the 
unification scale) when asymmetry, Ae N , is maximal. 

trivially extrapolated. The total number of events observed will be an interplay between 
the branching ratios and the production cross section. However, in the case of branching 
ratios, each of the decays, t\ —> xfb, t\ — > xtb and t\ —> Xit has a different reconstruction 
efficiency and asymmetry dilution that needs to be calculated. For example, we see from 
Fig. 5(b) that the branching ratio for the decay t± — > x^b increases noticeably as we vary 
4>A t from (ft A t — to 4>A t = I 71 "! due to this decay becoming kinematically more favourable. 
The xt has a large number of final states with no lepton however, so consequently the 
number of signal events decreases. Also, the xt decays generally contain extra jets that 
make the reconstruction of the event more difficult and thus reduce the efficiency of this 
channel. 

Figure 6(a) shows the asymmetry when all t± decay channels are considered and an 
estimate of the amount of luminosity required for a 3<r-observation (statistical errors 
only) of a non-zero asymmetry for pure tit\ production at the LHC We can see that the 
asymmetry is slightly diluted when all ti decay modes are included from I^^Jmax ~ 15% 
to l^l^lmax ~ 12.5%. The dilution is due to reconstructed events that are not originating 
from the signal process, Eq. (19). These events have no overall asymmetry and therefore 
simply dilute the signal. The horizontal lines show the estimate of the required luminosity 
required to see a certain asymmetry; an asymmetry can be seen at the 3<r level where 
the asymmetry curve in Fig. 6(a) lies outside the luminosity band. The luminosity bands 
are not flat because as discussed before, both the branching ratios and production cross 
section of the t\ vary with the phase 4>A t - We can see that in our scenario for pure t\t\ 
production, we expect a sensitivity for 0.57T < 0A t (mod it) < 0.9n with 500 fb _1 . With a 
combined analysis of both ATLAS and CMS data, this luminosity can be expected to be 
reached in the early 2020's if the LHC operates as is currently planned [69]. 

We can see the effect of varying the mSUGRA parameters tan/3 and Aq in Fig. 6(b). 
It is shown that as the value of either tan/3 or Aq is increased, we require more luminosity 
to see a statistically significant observation even with maximum asymmetry. An increase 
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Figure 7: Minimum luminosity required for 3 a- discovery in mo, m 1( / 2 plane (at unifica- 
tion scale) when asymmetry, Ae N , is maximal. Pure tit* production, all decay channels 
included, see Tab. 5 for branching ratios for the specific parameter point and Fig. 5 for 
how these alter with 4>A t - Purple area is ruled out by LEP direct detection [70] and red 
area has no two body decay X2 — * (a) With r tagging, (b) Without r tagging. 

in tan/3 decreases the sensitivity because the branching ratio x\ ~~ >* & is reduced. The 
reduction is due to f's becoming more mixed which increases the left handed component 
in the lighter f . Therefore, the f\ couples more strongly to the predominantly wino x\ 
and begins to dominate this decay channel at the expense of the signal process. A rise in 
A decreases sensitivity mainly because the CP-asymmetry is reduced. The reason is that 
after RGE running, an increase in Aq reduces the magnitude of the trilinear coupling A t 
that contains the phase, <pA t that we are interested in. Hence the CP effects are reduced. 

Similarly, Fig. 7(a) shows the effect of varying the mSUGRA parameters m and 
on the minimum luminosity required for an observation of CP effects. We note as a general 
trend that as m^j is increasing, we need more luminosity to observe the CP- violating 
triple products. This is due to the increase in t\ mass which reduces the production cross 
section for tit*. If we increase mo we see that a large area of the parameter space has no 
two body decay X2 ~~ ^ ± ^ T as m i± > m x°- 

Fig. 7(b) indicates the effect of having no hadronic r-tagging for the decay xt ~~ ^ ^\ v r- 
The r final state dominates the xt decay which in turn is the dominant product of the ti 
in low mass mSUGRA scenarios, Tab. 5. As stated in the beginning of Sec. 5 we assume a 
40% r-tagging efficiency and without this we lose approximately a factor of 2 in effective 
luminosity for our signal process. 

5.4 Impact of momentum reconstruction on SUSY background 
separation 

All of the previous section's results have assumed that the tit* process can be isolated 
effectively. However, in the mSUGRA scenarios investigated many other SUSY particles 
will be produced. Table 7 shows that the total production cross section for SUSY is ~ 25 
times greater than for titj production and we can therefore expect sizable backgrounds. 
We can also expect that the vast majority of the SUSY background processes will have 
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SUSY 


fit; Signal / SUSY Background 


Cross Section (pb r ) 


3.44 


80.1 




Events with 500 fbr 1 


1.7xl0 6 


4xl0 7 




Events with 500 fb" 1 
Initial selection 


32389 


410735 


0.079 


Events with 500 fb" 1 
Top Reconstruction 


7117 


64729 


0.11 


Events with 500 fb" 1 
Kinematic Reconstruction 


1213 


3759 


0.32 


Events with 500 fb" 1 
Extra SUSY cuts 


901 


967 


0.93 



Table 7: Cross section, number of events and signal to background ratio at the LHC with 
y/s = 14 TeV at LO for both the production channel t]i\ and inclusive SUSY production. 
All cross sections were calculated using Herwig++ [55,56]. 



no other spin correlated CP-sensitive triple product with the same final state and will 
therefore just act as a dilution to the CP-asymmetry by contributing to the denominator 
of Eq. (15). 

Table 7 shows that after the initial event selection and top reconstruction, the SUSY 
background is still ~ 10 times larger than the signal process. Note that if we apply 
the kinematical reconstruction to these events we see that we substantially reduce the 
background to be only rs 3 times larger. 

In order to observe CP-violating effects in tit* production at the LHC, however, the 
signal to background ratio may still be too high and consequently we need further cuts to 
isolate the signal process. We notice that in mSUGRA scenarios, the largest background 
comes from g production followed by the dominant decay to either sbottom, g — > b{b with 
a branching ratio of ~ 30%. The 6, decays dominantly to x^b or xt t which leads to a very 
similar final state as the signal process when combined with the opposite decay chain. 
The difference between the SUSY background and the tiS is that the g and first and 
second generation q have a higher mass. In addition, a gluino has in general one more 
decay vertex in the cascade decay producing another hard jet. These two factors mean 
that the average pt of the particles produced in the event will be higher and the number 
of jets will be greater, thus we can use these characteristics to discriminate the signal 
from the background. Hence we cut on the number of jets reconstructed in an event, 

Number of jets < 6. (33) 
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Figure 8: General SUSY production for the asymmetry Ae N . r tagging is included in 
both plots, (a) Asymmetry, Ae N , at reference point with 3a-luminosity lines shown, (b) 
Minimum luminosity required for 3<r-discovery in tan/3, A plane (at unification scale) 
when asymmetry, Ae N , is maximal. 



For the px cuts, we have, 

p T (Hardest Jet) < 200 GeV, (34) 

p T (2ndJet) < 130 GeV, (35) 

p T (3rdJet) < 80 GeV (if applicable), (36) 

p T (Any6Jet) < 150 GeV, (37) 

p T (Any Lepton) < 100 GeV. (38) 



Table 7 shows that after all these cuts are performed the signal to background ratio 
improves significantly and we now have roughly the same number of signal and background 
events in the sample. 

If we now re-evaluate the luminosity plots with the SUSY background included, 
Fig. 8,9, we see that more luminosity is now required to observe a statistically significant 
effect. Due to the background dilution of the asymmetry, we now have 1.4^1 max ~ 6.5% 
for our scenario Fig. 8. Consequently we are now only sensitive to phases between 
0.67r < (^(mod 7r) < 0.857T with 1 ab" 1 of data. If we look at the tan/3, A contour 
plot we see that sensitivity at the LHC for 1 ab -1 is only possible for small values of tan/3. 
A luminosity of 1 ab -1 would probably require an upgrade to the High-Luminosity LHC 
(HL-LHC) [71]. 

However, we would like to emphasise that it may be possible to substantially im- 
prove the statistical significance of an asymmetry measurement and return to close to the 
significance achieved when looking at a purely tit\ process, even with the same SUSY 
background. Namely, via measuring the SUSY spectra (in particular the g and b) a good 
estimate of the background should be possible. The background events can then be sub- 
tracted from the denominator of the asymmetry, Eq. (15), to give the true value of the 
asymmetry. Thus, the statistical significance should be much improved. 

We would also like to remind the reader that this subtraction only becomes reliable 
if the signal to background ratio is good enough, otherwise the signal is swamped by 
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Figure 9: General SUSY production for the asymmetry Ag N . Minimum luminosity re- 
quired for 3cr-discovery in m ,mi/ 2 plane (at unification scale) when asymmetry, Ae N , is 
maximal. Purple area is ruled out by LEP direct detection [70] and red area has no two 
body decay X2 ~~ >* ^ ± ^ T - (a) With r tagging, (b) Without r tagging. 



statistical fluctuations. Thus the momentum reconstruction procedure is vital since it 
significantly reduces the backgrounds that are present. 

Similarly, a more constrained area of observability is seen in the mo, mi/2 plane, 
Fig. 9(a). With 1 ab _1 of data, our study suggests that only if myi < 220 GeV will it 
be possible to observe a CP-phase in the stop sector. Again, we see the importance of 
r-tagging to our study from the difference between Fig. 9(a) and Fig. 9(b). If r-tagging 
is not used in the study, no CP- violation in the t\ sector can be observed with 1 ab _1 of 
data. 

5.5 Open experimental issues 

Although the presented study was completed at the hadronic level, a full detector simu- 
lation should be completed to confirm the conclusions of this paper. The most obvious 
experimental issue that could affect our results is the finite momentum resolution of the 
detector for both jets and leptons when performing momentum reconstruction. However, 
the resolution was tested with regards to momentum reconstruction in Ref. [29] with a 
significantly more complicated final state and it was found to have only a small effect. 

In terms of background suppression the mis-tagging of various objects could increase 
both the standard model and SUSY background. For the standard model background, 
the most obvious example is the it process generating a trilepton signal [63] . The process 
requires a jet to be mistagged as a lepton, which is not investigated in this study. The 
suitability of hadronic r-tagging in the study also needs to be investigated thoroughly 
as these are expected to have significant mis-identification rates [63]. However, this is 
beyond the scope of this theoretical study. 
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6 Conclusions 



In this paper we have investigated the problem of discovering CP-violating effects at the 
Large Hadron Collider. We studied t\t\ production and subsequent two body decays. 
Triple product correlations can be formed from the final state particles that are sensitive 
to the presence of complex phases in the model. Since triple products depend crucially 
on spin correlations and are therefore sensitive to CP-odd observables, they have been 
included both in the analytical calculation and the event generation, that has been per- 
formed using Herwig++ 2.3.2. The process of special interest in our case was the t\ decay 
into t and x° followed by two, 2-body leptonic decays. For this decay in our mSUGRA 
scenario one can expect an asymmetry in the triple product distribution of up to 15% 
when calculated in the rest frame of the decaying neutralino. The source of the CP viola- 
tion in our case was the phase of the trilinear coupling A t that attains a value of <f>^ t ~ 0.8 
when the asymmetry is maximum in our scenario. 

Due to the hadronic experimental environment of the LHC, precise measurements will 
be a challenge both from experimental and theoretical point of view. The rest frame 
CP-odd asymmetry is diluted by the high boosts of the produced particles and this makes 
an observation difficult. We studied the impact of momentum reconstruction of invisible 
LSPs to get access to the rest frame of the i\. Using a set of invariant kinematic conditions 
we showed that it is possible to fully reconstruct the production and decay process on an 
event-by-event basis. The reconstruction was performed on events including the parton 
shower and hadronisation. Having fully reconstructed events we are able to boost particle 
momenta back to the rest frame of the t\ and the maximum asymmetry is recovered 
to 15%. In addition, momentum reconstruction leads to a significant increase in the 
signal background ratio and thus is very important in attempting to isolate the process 
of interest. 

If we consider exclusive t\ production and all possible t\ decay chains the maximum 
asymmetry is diluted slightly to ~ 12.5%. In the mSUGRA scenario considered in this 
paper one should expect to see a 3a effect at C = 500 fb _1 for phases in the range 0.5 n < 
^(mod 7r) < 0.9 7r. If general SUSY production is considered, significant backgrounds 
to our signal process are present and extra kinematical cuts are required to remove this 
background. Even after these cuts some SUSY background remains and our maximum 
asymmetry is reduced to ~ 6.5%. To see a 3a effect at the LHC would require C = 1 ab _1 
of data for sensitivity to phases in the range 0.6 7r < 0^ t (mod n) < 0.85 ir. 

We emphasise that the asymmetry after momentum reconstruction is a much cleaner 
observable from a theoretical point of view, thanks to a well defined final state. Therefore, 
using the above technique provides prospects for the observation of CP-violating effects for 
a range of the phase <pA t after a few years of LHC running at the high luminosity. The full 
assessment of LHC's ability to resolve CP violation in the MSSM, however, will definitely 
require a detailed simulation of detector effects, SM and SUSY backgrounds which is 
beyond the scope of the present phenomenological analysis. The promising results of this 
study may encourage such further simulations. 
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Appendices 

A Mixing in the stop sector 

In the Minimal Supersymmetric Standard Model the stop sector is defined by the mass 
matrix M-i in the basis of gauge eigenstates (£l, t R ). The 2x2 mass matrix depends on 
the soft scalar masses Mq and My, the supersymmetric higgsino mass parameter /i, and 
the soft SUSY-breaking trilinear coupling A t . It is given as [72], 

M 2 (m\^m\ L m* LR m t \ 
1 \ m LR m t mf + m RR J 

where, 

m\ L = M| + m 2 z cos 2(3 (- - sin 2 6 W ) , (40) 
2 

m 2 RR = M? + -m| cos 2(3 sin 2 9 W , (41) 
m LR = A t - /I* cot (3 , (42) 

and tan/3 = v 2 /vi is the ratio of the vacuum expectation values of the two neutral Higgs 
fields which break the electroweak symmetry. From the above parameters only fi and A t 
can take complex values, 

A = |A|e^, /i=|/i|e^, (O<0^,0 M <2vr), (43) 

thus yielding CP violation in the stop sector. 

The hermitian matrix A^ 2 is diagonalized by a unitary matrix lZf, 

where we choose the convention m? < m? for the masses of t\ and to. The matrix 71+ 
rotates the gauge eigenstates, and tn, into the mass eigenstates t\ and t 2 as follows, 

ti\ ^ ( ' Il\ ( cos^f sin^ e -1 ^' \ f t L 



h)~ ni \t R J-\- sin 0- t &t cos 6- t ) { t R ) ' (45) 
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where 9f and fa are the mixing angle and the CP-violating phase of the stop sector, 
respectively. The masses are given by, 

m h,2 = \ ( 2m * + m lh + ™ 2 RR =F >J (rn? LL - m RR f + A\m LR \ 2 m^j , (46) 

whereas for the mixing angle and the CP phase we have, 

-m t \m LR \ 

cos 9- t = = , (47) 



nr. 



\m LR \ 2 + (m? - m\ L ) 2 



ti 
I 2 

sin 6~ t = , (48) 

Jm 2 \m LR \ 2 + (m| - m| L ) 2 

t ~ = arg(A-/i*cot/3). (49) 

By convention we take < 9f < n and < <pi < 2tt. It must be noted that <pi is an 
'effective' phase and does not directly correspond to the phase of any MSSM parameter. 
Instead, the phase will have contributions from both <pA t and </>^. However, in this study 
we set M = due to the EDM constraints. 

If m LL < m RR then cos 2 Q~ t > \ and t x has a predominantly left gauge character. On 
the other hand, if m LL > m RR then cos 2 6~ t < | and t\ has a predominantly right gauge 
character. 



B Mixing in the neutralino sector 



M 



N 



(50) 



In the MSSM, the four neutralinos X® (i — 1, 2, 3, 4) are mixtures of the neutral U (1) and 
SU(2) gauginos, B and W 3 , and the higgsinos, H® and H%. The neutralino mass matrix 
in the (B, W 3 , H%, H%) basis [73, 74], 

/Mi -m z Cf}S W m z Sf}S W \ 

M 2 mzcpcw -mzspcw 

—mzcpsw rrizCpCw — /i 

\ mzspsw -mzspcw —[i ) 

is built up by the fundamental SUSY parameters: the £7(1) and SU(2) gaugino masses 
Mi and M 2 , the higgsino mass parameter /i, and tan/3 = v 2 /vi (c@ = cos/3, sw = sin9 w 
etc.). In addition to the \i parameter, a non-trivial CP phase can also be attributed to 
the Mi parameter: 

Mi = \Mi\ e i0 \ (0 < 0i < 2tt) . (51) 

Since the complex matrix .Mat is symmetric, one unitary matrix N is sufficient to rotate 
the gauge eigenstate basis (B, W 3 , H 2 ) to the mass eigenstate basis of the Majorana 
fields x°i 

diag(m^o, m^o, m^o, m^o) = N*M.nN* , 



m xl < m x% < m 5& < m xl) 



"xi 



v x\> 



(52) 

The masses m^o (i = 1,2,3, 4) can be chosen to be real and positive by a suitable definition 
of the unitary matrix N. 
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C Interaction Lagrangian and couplings 

The interaction Lagrangian for the stop decay (tj — > x®t) is, 

A*x° = X°j KPl + kjPR) t 1* + h.c. , (53) 
where Pl,r — |(1 -F 7s)- The couplings are given by, 



, , -7= ^ f + ch/A^) - K t ^* 2 7V; 4 , (54) 



e _ f /l 

a. 



2\/2e 
3cw 



K = -i— n * N n ~ Y t n^Nj, , (55) 



where 1Z\j are the entries of stop mixing matrix, Eq. (45), and iV^ are the entries of the 
neutralino mixing matrix, Eq. (52). The top Yukawa coupling is given by, 

Y t = . (56) 

V2 mwsw sin (3 

The interaction Lagrangian for the neutralino decay (x° — > ££) is, 

£ «x° = 9fU I PR X? *l + gf Rj i Pl X°j Ir + h.c. (57) 

where g = ej %va.Q w . The couplings are given by, 

ft, = ^=(tan^7V. 1 + 7V. 2 ), (58) 

f R . = -V^tanewN*,, (59) 

(60) 

D Amplitude squared including full spin correlations 

D.l Neutralino production t\ — > x®t 

Here we give the analytic expression for the neutralino production density matrix [23], 
\M{h x°t)\ 2 = P(x° 3 t) + X P (x°j) + Z b P (t) + Efix^t) , (61) 

whose spin-independent contribution reads 

P(xjt) = (Kf + \bij\ 2 ){ptPy») - 2m t m^Re{a lj b\ j ) , (62) 

where pt and p^o denote the four- momenta of the t-quark and the neutralino xt- The 
coupling constants and bij are shown in Eq. (54,55) and by substituting the explicit 
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matrix elements of Eq. (45) we can show the specific parameter dependence [75], 



l a ijl + = 
= cos 2 9+ 



2 s 2 c 2 



+ 2 sin 9i cos 9i Yt 



-swNj! + c w N j2 + Y t 2 \N j4 \ z ) + sin 2 6- t 



Re 



q r 2 ' j 



V2 



SwCw 



e^~ -swN*, + c w N* 2 N j4 



N n \ 2 + Y 2 \N ]4 
2\/2e 



3c 



w 



Re[e- i<f >W n N; 4 ] 



(63) 



Re [aij&ij] = cos 2 9f —= Y t Re 

V2 " 



S\\rCw 



-swN^ + c w n; 2 iv; 4 



. 2a 2 ^ e 

+ sin 



* 3c 



YtReiN^N*, 



2 e 2 



+ sin 0? cos 0; Y? Re [e _i **"iV;|l =- Re 

1 L JJ 3s w c 2 w 



e^~ -swN*, + c w N* 2 N*, 



(64) 



The spin-dependent terms that depend on individual spin contributions are T-even 
and are given by, 



2 _ 






a ij 


2 _ 


a ij 





Z b P (t) = (\b^\ 2 -\a tJ \ 2 )m t (p^s b (t)), 



(65) 
(66) 



where (s b (t)) denote the spin-basis vectors of the neutralino x° (t-quark). Again 

the coupling constants can be expanded as, 



cos'% r/iiv^r - 



2 sin 0£ cos 9i Y t 



2 s 2 c 2 



+ sin 2 9- f 



8e 2 



q r 2 



V2 



Re 



SwCw 



+ 



N n \ 2 -Y 2 \N j4 
2y/2e 



3c 



Re^A^A^] 



(67) 



The terms that depend simultaneously on the spin of the top quark and of the neutralino 
can be split into T-even, ^ / p teven (x < jt) , and T-odd, ^ ( p t0 dd(x < jt) ■ The T-even contributions 
are as follows, 



E2Un(*?t) = 2Re{a l0 b: 3 )[{s\x^)Pt){s\t) P ^ - (m^O(^VW)] 



< * \ |7 „a/ ~0 



af~0\ n bf 



+ m t m^(s a (x {) j )s b (t))(\a ij \ 2 + \b 



■tj i 



(68) 



The T-odd contributions that generate the triple product correlations that we are inter- 
ested in are, 



Zf <o Ux J t) = -g 2 Im(a lJ b* j )f?, 
where the T-odd kinematical factor is given by, 

ft = w^(x>*°* 6 'W- 



(69) 
(70) 
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Sec. 2. 2 explains how this epsilon product generates the triple product observable. We 
again expand the coupling constant to see the functional dependence, 



Im [aij&*j] = cos 2 Qi 



SwCw 



Y t lm 



-swN*, + c w N* 2 N* 4 



+ sin 9~ t cos 9~ t (y? Im [e^N**] 
2 e 2 



3 SwC^y 



Im 



( -swN*, + c w N* 2 ) N;, 



(71) 



D.2 Neutralino decay — > ^ 

nalytical expressions for thi 

state £~ [76], 



We provide analytical expressions for the 2-body decay of the x° m t° a an< ^ the fi na l _ 



D(x 2) = g j\fL\ 2 Ho-ml}. 
The spin-dependent contribution is T-even and reads: 



D.3 Top decay t W + b 



(72) 



(73) 



We provide analytical expressions for the 2-body decay of the top quark into a FT-boson 
and the final-state bottom quark [77], 



D(t) = »-{m 2 t -2m 2 w + -f}. 
The spin-dependent contribution is T-even and reads: 



(74) 



Z» D (t) = -?-m t {(s b (t)p b ) + 



m 



At) Pw )}. 



(75) 
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